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Abstract 

Systems with long-range persistence and memory are shown to exhibit different precursory as well as 
recovery patterns in response to shocks of exogeneous versus endogeneous origins. By endogeneous, we 
envision either fluctuations resulting from an underlying chaotic dynamics or from a stochastic forcing 
origin which may be external or be an effective coarse-grained description of the microscopic fluctua- 
tions. In this scenario, endogeneous shocks result from a kind of constructive interference of accumu- 
lated fluctuations whose impacts survive longer than the large shocks themselves. As a consequence, 
the recovery after an endogeneous shock is in general slower at early times and can be at long times 
either slower or faster than after an exogeneous perturbation. This offers the tantalizing possibility of 
distinguishing between an endogeneous versus exogeneous cause of a given shock, even when there 
is no "smoking gun." This could help in investigating the exogeneous versus self-organized origins in 
problems such as the causes of major biological extinctions, of change of weather regimes and of the 
climate, in tracing the source of social upheaval and wars, and so on. Ref. p9| ] has already shown how 
this concept can be applied concretely to differentiate the effects on financial markets of the Sept. 11, 
2001 attack or of the coup against Gorbachev on Aug., 19, 1991 (exogeneous) from financial crashes 
such as Oct. 1987 (endogeneous). 



1 Introduction 

Most complex systems around us exhibit rare and sudden transitions, that occur over time intervals that 
are short compared to the characteristic time scales of their posterior evolution. Such extreme events ex- 
press more than anything else the underlying "forces" usually hidden by almost perfect balance and thus 

*We acknowledge helpful discussions and exchanges with Y. Malevergne, J.-F Muzy and V. Pisarenko. This work was partially 
supported by the James S. Mc Donnell Foundation 21st century scientist award/studying complex system. 
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provide the potential for a better scientific understanding of complex systems. These crises have funda- 
mental societal impacts and range from large natural catastrophes such as earthquakes, volcanic eruptions, 
hurricanes and tornadoes, landslides, avalanches, lightning strikes, meteorite/asteroid impacts, catastrophic 
events of environmental degradation, to the failure of engineering structures, crashes in the stock market, 
social unrest leading to large-scale strikes and upheaval, economic drawdowns on national and global scales, 
regional power blackouts, traffic gridlock, diseases and epidemics, and so on. It is essential to realize that 
the long-term behavior of these complex systems is often controlled in large part by these rare catastrophic 



events [25]. The outstanding scientific question is how such large-scale patterns of catastrophic nature might 
evolve from a series of interactions on the smallest and increasingly larger scales [p8|], or whether their origin 
should be searched from exogeneous sources. 

Starting with Hurst's study of 690 time series records of 75 geophysical phenomena, in particular river flow 



statistics, documenting the so-called "Hurst effect" of long term persistence [ ICQ, many studies in the last 
decades have investigated the existence of long memory effects in a large variety of systems, including me- 
teorology (wind velocity, moisture transfer in the atmosphere, precipitation), oceanography (for instance 
wave-height), plasma turbulence, solar activity, stratosphere chemistry, seismic activity, internet traffic, fi- 
nancial price volatility, cardiac activity, immune response, and so on. 

The question addressed here is whether the existence of long memory processes may lead to specific sig- 
natures in the precursory and in the relaxation/recovery/adaptation of a system after a large fluctuation of 
its activity, after a profound shock or even after a catastrophic event, that may allow one to distinguish an 
internal origin from an exogeneous source. Let us put this question in perspective with regards to the extinc- 
tion of biological species as documented in the fossil record. During the past 550 million years, there have 
been purportedly five global mass extinctions, each of which had a profound effect on life on Earth. The 
last end-Cretaceous mass extinction (with the disappearance of 39 — 47% of fossilizable genera and perhaps 
75% of species) marking the Cretaceous/Tertiary (K/T) boundary about 65 millions ago is often attributed 
to the impact of a huge meteor in the Yucatan peninsula [jl5|]. Another scenario is that a burst of active 
volcanism was the real origin of the extinction 0,|5|]. It has been suggested that this extinction was actually 
driven by longer-term climatic changes, as evidence by the fact that certain species in the Late Maastrichtian 
disappeared a distinctive time before the K/T boundary (JTJ, [T^]. A completely endogeneous origin has 
also been proposed, based on the concepts of nonlinear feedbacks between species [|2(], illustrated by 



self-organized criticality and punctuated equilibria ||2J, [24J] (see [ ]14| ] for a rebuttal). The situation is even 
murkier for the extinctions going further in the past, for which the smoking guns, if any, are not observable 
(see however the strong correlation between extinctions and volcanic traps presented in [||]). How can we 
distinguish between an exogeneous origin (meteorite, volcanism, abrupt climate change) and endogeneous 
dynamics, here defined as the progressive self-organizing response of the network of interacting species that 
may generate its own demise by nonlinear intermittent negative feedbacks or in response to the accumulation 
of slowly varying perturbations in the environment? Is it possible to distinguish two different exogeneous 
origins, one occurring over a very short time interval (meteorite) and the other extending over a long period 
of time (volcanism), based on the observations of the recovery and future evolution of diversity? 

The aviation industry provides another vivid illustration of the question on the endogeneous versus exoge- 
neous origin of a crisis. Recently, airlines became the prime industry victim of the September 11, 2001 
terrorist attacks. The impact of the downturn in air travel has been severe not just on the airlines but also on 
lessors and aircraft manufacturers. The unprecedented drop in air travel and airline performance prompted 
the US government to provide $5 billion in compensation and to make available $10 billion in loan guaran- 
tees. This seems a clear-cut case for an exogeneous shock. However, the industry was deteriorating before 
the shock of September 11. In the first eight months of 2001, passenger traffic for US carriers rose by an 
anemic 0.7 percent, a sharp fall from annual growth of nearly 4 percent over the previous decade [Qj, illus- 
trated by the record levels of the earned profits of $39 billion and of delivery of more than 4,700 jetliners 
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from 1995 to 2000. The US airlines' net profits dropped from margins of nearly 4 percent during 1998-2000 
to losses of greater than 3 percent during the first half of 2001, despite aggressive price cuts as airlines tried 
to fill seats and profits vanished. 

Many other examples are available. We propose to address this general question of exogeneous versus 
endogeneous origins of shocks by quantifying how the dynamics of the system may differ in its response to 
an exogeneous versus endogeneous shock. We start with a simple "mean field" model of the activity A(t) 
of a system at time t, viewed as the effective response to all past perturbations embodied by some noisy 
function r?(r), 



A(t) = f dr 7](t) K(t - t) 

J — oo 



(1) 



where K(t — r) can be called the memory kernel, propagator, Green function, or response function of the 
system at a later time t to a perturbation r/(r) that occurred at an earlier time r. Notwithstanding the linear 
structure of ([]]), we do not restrict our description to linear systems but take ([]]) as an effective coarse- 



grained description of possible complex nonlinear dynamics. For instance, it has been shown [30] that 



the extremal nonlinear dynamics of the Bak and Sneppen model and of the Sneppen model of extremal 



evolution of species, which exhibit a certain class of self-organized critical behavior [|27[], can be accurately 
characterized by the stochastic process called "Linear fractional stable motion," which has exactly the form 
(|l|) for the activity dynamics. 

Expression ([I]) contains for instance the fractional Brownian motion (fBm) model introduced by Mandelbrot 
and Van Ness [jl6|] as a simple extension of the memoryless random walk to account for the Hurst effect. 
From an initial value Bh(0), we recall that the fBm is defined by 

B " (t > - B »<°> = mwm L dT ,<T> K{t ~ T) ■ <2) 

where dr t](t) = dW T is usually taken as the increment of the standard random walk with white noise 
spectrum and Gaussian distribution with variance E[dW T ] = dr and the memory kernel Kit — r) is given 
by 

K(t-r) = (t-T) H ~i, forO <T<t (3) 
= (t-r) H "2 - (-t) h ~^ , forr<0. (4) 

For H > 1/2, the fBm Bh (t) exhibits long term persistence and memory, since the effect of past innovations 
of dW T is felt in the future with a slowly decaying power law weight Kit — r). For our purpose, the fBm is 
non-stationary and it is more relevant to consider globally statistically stationary processes. 

Here, we consider processes which can be described by an integral equation of the form ([]]) and (|2|) but 
with possibly different forms for the noise innovations r] and for the memory kernel K. Simple viscous 
systems correspond to K (t — r) oc exp[— (t — r) /T], where T is a characteristic relaxation time. Complex 
fluids, glasses, porous media, semiconductors, and so on, are characterized by a memory kernel K(t — r) oc 
g-a(t-T)' 3 ^ with < P < 1, a law known under the name Kohlrausch-Williams-Watts law jj23"l]. It is also 
interesting to consider fractional noise motion (fNm) defined as the time derivative of Bu(t), which does 
possess the property of statistical stationarity. A fNm is defined by (|l|) with 

Km ^- r)= {t _ r) l-H = (t- r y-°> (5) 

for H = 1/2 + 8. Persistence 1/2 < H < 1 (respectively antipersistence < H < 1/2) corresponds to 
< 8 < 1/2 (respectively — 1/2 < 8 < 0). Such a memory kernel describes also the renormalized Omori's 



law for earthquake aftershocks [26, 
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2 Exogeneous versus endogeneous shock 



In the following, we consider systems described by a long memory integral ([]]) with kernel K(t) decaying 
faster than at large times, so as to ensure the condition of statistical stationarity. This excludes the 
fBm which are non-stationary processes but includes the fNm. 

2.1 Exogeneous shock 

An external shock occurring at t = can be modeled in this framework by an innovation which takes the 
form of a jump Aq 5(t). The response of the system for t > is then 

A(t) = A K(t)+ f drr](r)K(t-r). (6) 

J — oo 

The expectation of the response to an exogeneous shock is thus 

E exo [A(t)]=A K(t) + n(r ) ), (7) 

where (77) is the average noise level and n = f Q +oc dr K(t) is the average impact of a perturbation which 
is usually smaller than 1 to ensure stationarity (this corresponds to the sub-critical regime of branching 
processes [Q]). 

The time evolution of the system after the shock is thus the sum of the process it would have followed in 
absence of shock and of the kernel K(t). The response Aq K(t) to the jump Aq 5{t) examplifies that K(t) 
is the Green function or propagator of the coarse-grained equation of motions of the system. Expression (^) 
simply expresses that the recovery of the system to an external shock is entirely controled by its relaxation 
kernel. 



2.2 Endogeneous shock 

2.2.1 Conditional response function 

Let us consider the natural evolution of the system, without any large external shock, which nevertheless 
exhibits a large burst A(t = 0) = Aq at t = 0. From the definition ([I]), it is clear that a large "endogeneous" 
shock requires a special set of realization of the innovations {rj(t)}. To quantify the response in such case, 
we recall a standard result of stochastic processes with finite variance and covariance that the expectation of 



some process X(t) conditioned on some variable Y taking a specific value Yq is given by JTT[ ] 



E[X(t)\Y = Aq] - E[X(t)} = (A - E[Y}) C ° V ^ 2 j' y) , (8) 

where E[Y 2 ] denotes the expectation of Y 2 , Cov(X(t),Y) is the covariance of X and Y, E[X(t)\ and 
E[Y] are the (unconditional) average of X(t) and of Y. Expression (||) recovers the obvious result that 
E[.X~(f)|Y = ^o] = E[X(t)] if X and Y are uncorrected. A result generalizing <jsj> holds when t/(t) has an 
infinite variance corresponding to a distribution with a power law tail with exponent smaller than 2 [gj. 

Let us assume that the process A(t) and the innovations r/'s have been defined with zero mean, which is 
always possible without loss of generality by a translation. Let us call X(t > 0) = A(t) and Y = ^4(0). 
Under the assumption that the noise t/(t) has a finite variance, we obtain from ([[]) 

Cov(A(t), A(0)) = ! dr K{t - r) K(-t) , (9) 
J —00 
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and 



E[A(0y] = / dr [K(-t)Y . (10) 



For stationary processes such that K(t) decays faster than 1/ y t so as to make the integral in (|10J) convergent, 
E[A(0) 2 ] is a constant. We thus obtain the posterior (t > 0) response (above the stationary average) to an 
endogenous shock occurring at time t = under the form of a conditional expectation of A(t), conditioned 
by the existence of this shock: 

r+oo 

E eado [A(t)\A(0) = A ] ocA du K{t + u) K{u) , (11) 

Jo 

for large Aq. This relaxation of the activity after an endogeneous shock is in general significantly different 
from that given by (^) following an exogeneous shock. 

2.2.2 Conditional noise trajectory 



What is the source of endogeneous shocks characterized by the response function Qllp? To answer, let us 
consider the process W(t) = f_ dr t)(t), where fj(t) = r)(t) — (rj) defines the centered innovations 
forcing the system ([!]). Using the property t^>, we find that for t < 

E endo [W(t)\A(0) = A ] = CW ^ r ( * ) ; n f 1 ( ° )] • (4> - E[A\) oc (A - E[A}) f dr K(-r) , (12) 

where E endo [W(t) \A(0) = A ] = for t > since the conditioning does not act after the shock. Expression 
(|l2|) predicts that the expected path of the continuous innovation flow prior to the endogeneous shock (i.e., 
for t < 0) grows like AW(t) = fj(t)At ~ K(—t)At upon the approach to the time t = of the large 
endogeneous shock. In other words, conditioned on the observation of a large endogeneous shock, there is 
specific set of trajectories of the innovation flow rj(t) that led to it. These conditional innovation flows have 
an expectation given by (|l2|). 

Inserting the expression ( |l2| ) for the average conditional noise in the definition ([[]) of the process, we obtain 
an expression proportional to <JTT|). This shows that the precursory activity preceding and announcing the 
endogeneous shock follows the same time dependence as the relaxation ( |TT| ) following the shock, with the 
only modification that t (for t > counting time after the shock at t = 0) is changed into —t (for t < 
counting time before the shock at t = 0). We can also use (12) into ([I]) and calculate the activity after the 



endogeneous shock to recover (11). These are two equivalent ways of arriving at the same result, the one 



using ( |12[ ) illuminating the fundamental physical origin of the endogeneous response. 



These results allow us to understand the distinctive features of an endogeneous shock compared to an exter- 
nal shock. The later is a single very strong instantaneous perturbation that is sufficient in itself to move the 
system significantly according to (|6|). In contrast, an "endogeneous" shock is the result of the cumulative ef- 
fect of many small perturbations, each one looking relatively benign taken alone but, when taken all together 
collectively along the full path of innovations, can add up coherently due to the long-range memory of the 
dynamical process. In summing, the term "endogeneous" is used here to refer to the sum of the contribution 
of many "small" innovations adding up according to a specific most probable trajectory, as opposed to the 
effect of a single massive external perturbation. 

2.3 Numerical simulation of an epidemic branching process with long-range memory 

To illustrate our predictions (0) and (|TI|), we use a simple epidemic branching model defined as follows. The 
model describes the time evolution of the rate of occurrence of events as a function of all past history. What 
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is called "event" can be the creation of a new species or a new family of organisms as in [m, the occurrence 
of an earthquake as in [|^, §], the amplitude of the so-called financial volatility as in [ |29| ] or of aviation 
traffic, a change of weather regime, a climate shift and so on. The rate \(t) of events at time t is assumed to 
be a function of all past events according to 

i | ti<t 

where the sum is carried over all past events that occurred at times t{ prior to the present t. The influence 
of such an event at a previous time ti is felt at time t through the bare propagator <j)(t — U). In our present 
illustration, we consider a process equivalent to a fNm with Hurst exponent H = 1/2 + 9, which can be 
shown to correspond to the choice 4>{t) = 9 c e /(t + c) 1+e , where c is an ultra-violet regularization time 
embodying a delay process at early times in the activity response after an event. Indeed, the Master equation 



corresponding to the process ( |13| ) can be shown [ ]26| , to be nothing but (|1|) with the renormalized or dressed 
propagator K(t) oc l/(f + c) l ~ e . 

Numerical simulations of the epidemic branching process are performed by drawing events in succession 
according to anon-stationary Poisson process with instantaneous rate (13). Figures [l] and ^] show successive 



magnifications of time series of the activity rate after an exogeneous shock and around an endogeneous 
shock, respectively, in order to visualize the precursory and relaxation activities. In figure ||, an external 
source of activity necessary for seeding has been added as a Poisson process of rate fi = 10~ 3 corresponding 
on average to one external event over a time interval of 1000. The most striking visual difference is the 
existence of the precursory signal occurring at many time scales for the endogeneous shock. Figure || 
quantifies the precursory and relaxation rates associated with activity shocks. The top panel shows the 
relaxation of the activity (rate of events) following an external shock compared to that after an endogeneous 
shock, for a single realization. t c is the time of the shock. The horizontal axis is t — t c for the relaxation 
of the activity after the shock. The precursory activity prior to the shock is also shown for the endogeneous 
shock as a function of t c — t. The bottom panel shows the same three activity functions after averaging 
over many realizations, translating time in the averaging so that all shocks occur at the same time denoted 
t c . The prediction (0) states that the relaxation of the activity after an exogeneous shock should decay as 
K(t) oc l/(t — t c + c) 1-0 while the decay after an endogeneous shock should be given by ( |I~ij ) which 
predicts the law oc l/(t — t c + c) 1-26 *, that is, a significantly smaller exponent for 9 > 0. Similarly, we 
predict that the precursory activity prior to an endogeneous shock should increase as oc l/(t c — t + c) 1 ^ 29 . 
These predictions are verified with very good accuracy, as seen in figure |3} 

These simulations confirm that there is a distinctive difference in the relaxation after an endogeneous shock 
compared to an exogeneous shock, if the memory kernel is sufficiently long-ranged. For a single realization, 
there are unavoidable fluctuations that may blur out this difference. However, we see a quite visible precur- 
sory signal (foreshock activity) that is symmetric to that relaxation process in the case of an endogeneous 
shock. This follows from the model used here which obeys the time-reversal symmetry. This may be used 
as a distinguishing signature of an endogeneous shock. 



3 Classification of the distinctive responses for different classes of memory 
kernels 

The family of power law kernels used in the simulations presented in figures [|, H and [3] are only one pos- 
sibility among many. Our formalism allows us to classify the distinctive properties of the relaxation and 
precursory behaviors that can be expected for an arbitrary memory kernel. We now provide this classifica- 
tion by studying (0) and (11). 
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3.1 Short-time response 



We compare the initial slopes of the relaxations after the occurrence of the shock at t = 0. Thus, by 
short-time, we mean the asymptotic decay law just after the shock. For this, we expand ([7|) to get 



E m [A(t)] = A K(0) 



A K(0) 



dhxK, 
1 H — U=o 



dt 



t + 0(t 2 



where K'{t) denotes the derivative of K(t) with respect to time. 
Similarly, expanding the integral in dill ) for short times, we obtain 

ldlnF, 



E endo [A(t)\A(0) = A ] cx A o F(0) 



1 + 



2 dt 



t + 0(t 2 



where 



/+oo 
du[K(u)f 

is a monotonically decreasing function of time. 
It is convenient to use the parameterization 



F(t) 



-g(t) 



where g(t) is an monotonously increasing function of time. Inserting ( fl7| ) in ([T^) and (15) leads to 



E cxo [A(t)] = A Q K(0) 



and 



E endo [A(t) \A(Q) =Aq] oc AoF(0) 



l-lg'(0)t + O(t 2 ) 



(14) 



(15) 



(16) 



(17) 



(18) 



(19) 



1. For g"(0) = 0, that is, git) = 2at corresponding to an pure exponential relaxation Kit) oc exp[— at], 
the velocities of the responses to an exogeneous and to endogeneous shock are identical; 

2. for g"(0) > corresponding to a super-exponential relaxation K it) oc exp[— at c ] with c > 1, the 
exogeneous relaxation is slower than the endogeneous one; 

3. for g"(0) < corresponding to a sub-exponential relaxation such as a stretched exponential Kit) oc 
exp[— at c ] with c < 1 or to the family of regularly varying functions such as power laws, the exoge- 
neous relaxation is faster than the endogeneous one. 



The exponential relaxation thus marks the boundary between two opposite regimes. As is intuitive, a sub- 
exponential relaxation betraying a long memory process leads to a slower short-time recovery after an en- 
dogeneous shock, because it results from a long preparation process (|l2|). 



3.2 Asymptotic long-time response 

Since K(t) is a monotonously decaying function, Kit + u) < K(t) for any u > 0. This leads to the 
following inequality 

r+oa 

E endo [^(t)|A(0) = Aq] < A Q K{t) / du K(u) , (20) 

J o 
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which is valid if the integral Jq + °° du K[u) exists, that is, if K (t) decays faster than 1/t at large times. This 
shows that, as soon as K(t) -C C/t for any positive constant C, E ea do[^(i)|^(0) = A)] < E cxo [A(t)]. 
But the difference may be small and inobservable. For instance, for K(t) oc l/t 1+e with 9 > 0, a careful 



examination of the integral in ( Jl 1| > shows that, due to the contribution of the conditional noise close to the 
shock, we have 

E endo [A(t)\A(0) = ,4 ] oc ~ E cxo [A(t)} . (21) 

Thus, there is no qualitative difference in the relaxation rates of an endogeneous shock and exogeneous shock 
in this case: the contributions of all the conditional activity prior to the endogeneous shock is equivalent to 
that of the shock itself. A more elaborate and analysis specific to the problem at hand must be performed to 
predict the prefactors that will be different in the endogeneous and exogeneous cases. 

In constrast, for memory kernels K(t) oc 1/t 1 ^ 6 with 9 > decaying slower than 1/t, as for a stationary 
fNm of the form ©, we obtain 

E cndQ [A{t)\A{$) = A ] oc > E exo [A(t)} oc . (22) 

In this case, the relaxation following an endogeneous shock decays significantly more slowly than for an 
exogeneous shock. This case is examplified in figure ||. In the long time limit, the decay law 1/t thus marks 
the boundary between two opposite regimes. 



3.3 Illustration 

An illustration of this critical behavior is provided by the response of the price volatility a At at scale At 
defined as the amplitude (absolute value) of the return rjjt(t) = ln\p(t)/p(t — At)] = e(t) ■ 0At(i) = 
e(t) ■ e WA '(*). of a financial asset. e(t) is a random sign. Indeed, financial price time series have been shown 



to exhibit a long-range correlation of their log-volatility a; At. described by a model Q22| , |29| ] in which uJAt(t) 
follows the process ([]]) with 



K At (t) ~ K J — for At < t < T , (23) 

where T w 1 year is a so-called integral time scale. This form (^3|) corresponds to the parameterization (§) 
with 9 = 1/2. Sornette et al. [ [2^ ] have shown that there is a clear distinction between the relaxation of stock 
market volatility after an exogeneous event such at the September 11, 2001 attack or the Aug., 19, 1991 
coup against Gorbachev and that after an endogeneous event such as the October 19, 1987 crash. In this 
model, the long-range memory acting on the logarithm of the volatility induces an additional effect, namely 
the exponent of the power law relaxation after an endogeneous shock is a linear function of the amplitude 
of the shock. 



3.4 Synthesis of the asymptotic short- and long-time regimes 

We have found two special functional forms for the response kernel K(t) oc 1/t and K(t) oc exp[— at], 
which are "invariant" or indifferent with respect to the endogeneous versus exogeneous origin of a shock. 
Thus, for normal exponential relaxation processes as well as for power relaxation oc 1/t, the functional form 
of the recovery does not allow one to distinguish between an endogeneous and an exogeneous shock. 

These two invariants K(t) oc exp[— at] and K(t) oc 1/t delineate two opposite regimes, the first one for 
short-time scales and the second one for long-time scales: 
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1. for K(t) oc l/t 1 with 9 > 0, the endogeneous response decays more slowly than the exogeneous 
response, at all time scales; 

2. for exp[— at] <C K(t) <C l/t for any positive a, the endogeneous response decays more slowly than 
the exogeneous response at short time scales and has the same dependence as the exogeneous response 
at long time scales; this regime describes for instance the stretched exponential relaxation of complex 
fluids alluded to above; 

3. for K(t) <C exp[— at] for any positive a, the endogeneous response decay faster than the exogeneous 
response at all time scales. 



More complicated behaviors can occur when the memory kernel K(t) exhibits a change of regime, crossing 
the exponential and/or l/t boundaries at certain time scales. Each situation requires a specific analysis 
which yields sometimes surprising non-intuitive results Qj. 



4 Conclusion 



We think that the conceptual framework presented here may be applied to a large variety of situations, 
beyond those alluded to in the introduction. For instance, the result ( fl2| ) has been shown to explain the 
so-called inverse Omori's law for earthquake foreshock activity before a mainshock, in a simple model of 
earthquake triggering [g]. The same mechanism may explain the premonitory seismicity pattern known as 
"burst of aftershocks" [O]: a mainshock with an abnormally large number of aftershocks has been found to 



be a statistically significant precursor to strong earthquakes [21]. 



Many dynamical systems in Nature, such as geophysical and biological systems (immune network, memory 
processes in the brain, etc.), or created by man such as social structures and networks (Internet), States 
and so on, exhibit long-memory effects due to a wealth of possible mechanisms. For instance, Krishan 
Khurana at UCLA has suggested to us that the concept proposed here could explain that endogeneous civil 
wars have long-lasting effects with slow reconstruction compared with the fast recovery after exogeneous 
wars (that is, imposed or coming from the outside). The increasing emphasis on the concepts of emergence 
and complexity has emphasized an endogeneous origin of the complicated dynamical behavior of complex 
systems. In reality, most (so-called) complex systems are the result of their internal dynamics/adaptation 
in response to a flow of external perturbations, but some of these external perturbations are rare extreme 
shocks. What is the role of these exogeneous shocks in the self-organization of a complex system? Can 
one distinguish the impact of extreme exogeneous shocks from an endogeneous organization at different 
time scales? Our present analysis has just scratched the surface of these important and deep questions by 
suggesting an angle of attack based on the conditional historical process at the basis of strong endogeneous 
fluctuations. Extensions of the present simplified framework involve the generalization to multidimensional 



coupled processes such as in [12] and to nonlinear spatio-temporal processes. 



We are grateful to A.B. Davis. V. Keilis-Borok and V.F Pisarenko for useful exchanges. This work was 
partially supported by the James S. Mc Donnell Foundation 21st century scientist award/studying complex 
system. 
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Figure 1 : Rate of activity following an exogeneous shock in a numerical simulation of the epidemic branch- 
ing model (13) generated with a memory kernel decaying as a power law <ft(t) ~ l/(i + c) 1+e with pa- 
rameters 9 = 0.1 and c = 0.001 without a constant source term (fj, = 0). The rate of activity following 
an exogeneous shock that occurred at t = is shown at increasing magnification from top to bottom. It is 
evaluated using a bin size decreasing by factors of 10 from St = 10 (a) to St = 0.01 (d). Averaging over 
many such realizations would yield the average power law decay K(t) ~ l/(i + c) 1-61 predicted by (f7|). 
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Figure 2: Rate of activity prior to and after an endogeneously generated major burst of activity generated 
by a numerical simulation of the epidemic branching model (13) with a power law kernel with the same 
parameters 6 = 0.1 and c = 0.001 as in figure [l] with in addition a constant Poisson source term with rate 
fi = 0.001 corresponding, on average, to one event added from an external source per 1000 time units. Most 
of the observed activity is thus the result of interactions between events. The rate of activity close to the 
largest peak of activity is shown at increasing magnifications from top to bottom and is evaluated as in figure 
[j]. Both precursory and relaxational processes can be observed at many time scales. 
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Figure 3: Rate of activity for a single sequence (a) of the epidemic branching model defined by ( |l3"| ) gener- 
ated with a memory kernel decaying as a power law of time with the parameters 9 = 0.2 and c = 1CT 3 and 
averaged over many sequences (b). The exogeneous relaxation is shown with diamonds, the endogeneous 
relaxation is shown as crosses and circles are for the precursory activity in the endogeneous case. Large 
fluctuations are observed in the precursory activity and in the endogeneous relaxation when looking at a 
single sequence, due to the small number « 100 of observed events. Averaging over 50 realizations, we 
see clearly the faster decay rate ~ \/t l ~ e for the exogeneous relaxation predicted by (0) compared with the 
endogenous one predicted by (fTT|). The same decay rate ~ l/i 1-26 * predicted by (jXTJ) is observed for both 
the endogeneous precursory and post-event relaxation. 
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